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A class of inflation theories called a-attractors has been investigated recently with interesting 
properties interpolating between quadratic potentials, the Starobinsky model, and an attractor 
limit. Here we examine their use for late time cosmic acceleration. We generalize the class and 
demonstrate how it can interpolate between thawing and freezing dark energy, and reduce the fine 
tuning of initial conditions, allowing w « —1 for a prolonged period or as a de Sitter attractor. 


I. INTRODUCTION 


Scalar fields play a fundamental role in cosmology, with 
the Higgs field in the standard model of particle physics, 
the inflaton seeding early universe quantum perturba¬ 
tions that develop into cosmic structure and a primor¬ 
dial gravitational wave background, and possibly a dark 
energy field that is responsible for late time cosmic accel¬ 
eration. One of the exciting recent developments in in¬ 
flation is the discovery of classes of theories that connect 
different scalar field theories. Here we examine aspects 
of this in the context of late time dark energy. 

The inflationary developments involve conformal ^- 
attractors and Kahler curvature a-attractors. These have 
the property of interpolating between power law scalar 
field potentials or, e.g., Higgs inflation, Starobinsky in¬ 
flation, and a fixed point [1-3]. Thus such theories make 
definite predictions for the primordial scalar perturba¬ 
tion tilt Ug and the tensor (gravitational wave) to scalar 
power ratio r, both measurable by cosmic microwave 
background experiments, with all models approaching a 
particular limit. 

We focus here on the minimally coupled a-attractors 
and investigate the dark energy properties of such mod¬ 
els. Several of the forms within this class have high en¬ 
ergy physics motivations, from supergravity and other 
theories. It would be interesting to evaluate whether 
these models have characteristics useful for late time ac¬ 
celeration, and possibly an improvement over standard 
quintessence potentials. 

In Sec. H we lay out the basic of the models and gen¬ 
eralize them to a family of models. Section HI addresses 
the attractor behavior and presents numerical results for 
the dark energy equation of state evolution. We compare 
generalized a-models to standard quintessence in Sec. IV 
and discuss their advantages, summarizing and mention¬ 
ing future work in Sec. V. 


II. a-MODELS 
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where Mp is the Planck mass, a is a parameter and a/^ 
is the potential function. The Starobinsky model corre¬ 
sponds to a = 1. 

Through a held redefinition 


(/> = \/6a tanh , (2) 


the kinetic term becomes canonical and the potential 
function 


= af . (3) 

Note that the function / is not wholly arbitrary since the 
field redefinition breaks down at ^ = oo. 

For example, if we tried to make a standard 
quintessence model be described by an a-model, this is 
not generally possible. Consider dark energy with a con¬ 
stant equation of state (pressure to energy density) ratio 
w. This is given by a potential [4] 

V ^ sinh“^^^'''“^^'’"'(^ — ^*) [constant w] . (4) 

and is accommodated by a function 

-{l+w)/\w\ 
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where x = tanh {cj)/ v^fe). However this function vanishes 
at (j) = oo (a: = 1), and furthermore blows up at <)> = 0. 

For inflation, two functional forms have been used, the 
T-mo del 



f{x) = cx , (6) 

and the Starobinsky form (with the Starobinsky model 
having a = 1) 
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The a-model as written for inflation has a non- 
canonical kinetic term and a potential, with Lagrangian 
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Note that they are equivalent at small x, and because 
X (X </) at small (j) then near the origin both models give 
quadratic potentials. If we were to take f{x) oc for 
X <C 1 then we could match onto any monomial potential 
(see, e.g., [3]). However here we will instead generalize 
to interpolating and extrapolating these two models, to 
study the deviations from quadratic behavior (but see 
Sec. IV). Our generalized a-model takes 


f{x) 
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The quantity c scales the amplitude of the potential, and 
can be fixed by requiring a certain dark energy density to¬ 
day; that is, c is effectively equivalent to the present 
dark energy density in units of the critical density. The 
parameter a scales the field value (f). 

As X —> 1, its maximum value, the potential goes 
to a constant. Basically it plateaus at a level Vp = 
ac^2“^". Thus the generalized a-model potential looks 
like a quadratic potential at small </), out to a width 
(j) ~ ^/a, and then flattens to a plateau at large 4>. Fig¬ 
ure 1 shows this potential for values of n = 0-3, and the 
comparison to a standard quadratic potential at small (j). 
(The potential is symmetric about x = 0 but we show 
only the positive half since the dynamics is symmetric.) 
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FIG. 1. The potential is plotted as a function of the field 4> in 
units of ^/a, for various values n. At small field values 4> ^ 
ydJ the a-model resembles the quadratic potential V oc 
(short, dotted green curve), while at large values it flattens 
to an uplifted plateau. 


III. EVOLUTION AND ATTRACTORS 


In the a-attractor family of inflation theories, the word 
“attractor” refers to their common values for the scalar 
perturbation power spectrum tilt Us and tensor to scalar 
power ratio r in the limit a —> 0, while for large a the 
models resemble monomial potentials 4>p [2]. Here we 
are interested in the dynamics of the scalar field, and its 
equation of state evolution, in particular whether it gives 
rise to late time cosmic acceleration. 

From Fig. 1 we can already expect the answer: if the 
field rolls into the quadratic region of the potential, then 
it should act like dark energy from a quadratic potential 
at late times, while if it remains on the plateau then it 
should act like a cosmological constant. Looking at this 
a bit more carefully, we see that the potential for large (j) 
is 
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This is basically an uplifted exponential potential. The 
exponential potential was one of the original dark energy 
models [5]. Moreover, the quadratic potential is in the 
thawing class of dark energy, moving away from a cosmo¬ 
logical constant state, while the exponential potential is 
in the freezing class, moving toward a cosmological con¬ 
stant. (See [6, 7] for discussion of these two main classes 
of dark energy.) 

Thus these a-models appear to conjoin in a way these 
two classes. Furthermore, because the transition depends 
on the value of a, we see that large a moves the plateau 
(and hence freezing behavior) to larger and larger </), es¬ 
sentially shrinking the plateau to a point maximum and 
making the potential look like a hilltop potential [8]. 

In fact, while these behaviors effectively hold, there is 
a formal de Sitter fixed point at w = —1 only for a par¬ 
ticular range of n (though for other n the field has nearly 
this behavior, for tens of e-folds or more, as we discuss 
below). Therefore we use the terminology a-model rather 
than a-attractor, in the dynamical sense. However, in the 
inflationary sense of a common limit for large a models, 
and for small a models, there are attractors as we show 
in Sec. IV. 

We numerically solved the scalar field equations of mo¬ 
tion and Friedmann equations for the background expan¬ 
sion to verify these behaviors. Indeed, for initial field val¬ 
ues (j)i the field gradually thaws and rolls down the 

potential, growing from equation of state ratio rc = — 1 to 
less negative values. It eventually (perhaps in the far fu¬ 
ture) oscillates around the quadratic minimum, giving a 
time averaged equation of state ratio {w) = 0. However, 
it can produce a long (if temporary) period of cosmic ac¬ 
celeration. Figure 2 shows the present equation of state 
ratio wq as a function of (j)i for various values of n. 

For (j) ^ i/a, one might suspect that since the plateau 
is not perfectly flat, but tilted upward, that the field 
could slide down into the minimum. However, since the 
field kinetic energy is rapidly damped away in an accel¬ 
erating universe, this is essentially negligible. Indeed, 
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FIG. 2. The present value of the equation of state ratio wo 
is plotted as a function of the initial field value in units 
of ^/a, for various values of n (solid curves). For contrast, 
this relation is also shown for a quadratic potential (dotted 
green curve, where the x-axis is simply For a given value 
4>i, the a-model can achieve wq closer to —1 (as preferred by 
observations). Conversely for a given bound on the deviation 
of Wo from —1, the a-model can be less fine tuned in the initial 
field value. 


traditional skating fields [7] where V = 0 have kinetic 
energy vanishing as a“®, where a is the cosmic expan¬ 
sion factor; fields on the plateau could be called elevated 
skaters. For example, for n = 0, where Fig. 1 shows the 
plateau begins at (^ > 5i/a, the equation of state evolves 
from Wo = —0.9995 today to only w(a = 7.5) = —0.9988 
for (j)i = 5^/a. If (pi = 20^/a, then the field rolls only 

AcP/Mp = 6 X 10-^. 

(One might imagine stranger situations, in which the 
field starts near the minimum but with initial velocity 
away from it high enough to reach the plateau, or starting 
on the plateau and shooting through the minimum to 
reach the x < 0 plateau. However, these are basically 
moot because the rapid Hubble expansion at early times 
redshifts away any large kinetic energy, requiring extreme 
fine tuning for these situations to occur.) 

Our generalization shows further interesting proper¬ 
ties. Figure 1 illustrates that the plateau begins at 
smaller p as n increases. Indeed, for n = 2 and pi = 3\/a, 
the equation of state evolves from wq = —0.9999 to 
w{a = 7.5) = —0.9998. Notably, for n > 2 we have a 
true attractor since the potential has a maximum on the 
plateau (at x = l/(n— 1)), and slopes down, rather than 
up, to its asymptotic value as seen from Eq. (9). For ex¬ 
ample, for n = 3 and pi > 1.35-^0 the asymptotic future 
state is de Sitter, w = — 1. 


The potential slope V = dV/dp is given by 
2 - x)[l + {1 - n)x] 

Vs (l + x)2" ■ ^ ’ 

Figure 3 shows the slope for various values of n. We see 
that the slope rapidly approaches zero for large p, mak¬ 
ing the plateau quite flat. Recall from the Klein-Gordon 
equation of motion for the scalar field that once the slope 
is negligible then the field quickly freezes in place. For 
example, for n = 1 and pi = 5-^0, the field rolls less 
than Ap/Mp = 0.06 in its entire history. Thus even 
when n < 2 and there is no formal de Sitter attractor, 
one can still have w r; — 1 for a long time. 



FIG. 3. The first derivative of the potential is plotted as a 
function of the field p in units of y/a, for various values n 
(increasing from top to bottom curves). Note that all deriva¬ 
tives approach zero at large field values (since the potential 
plateau is asymptotically flat). 


IV. RELATION TO QUINTESSENCE MODELS 

At small p, the generalized a-model reduces to the 
quadratic potential, 

v{p^V^)^^p^. ( 11 ) 

oa 

That is a thawing dark energy model, and to keep the 
present equation of state ratio wg near enough to —1 to 
be compatible with observational constraints, the initial 
position of the field requires some fine tuning. In partic¬ 
ular, starting the field too close to the minimum means 
that it will never achieve a present dark energy density in 
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units of the critical density = 0.7. For example, in 
order to reach this requires </), > 1.42 Mp, while to also 
keep wq < —0.9 requires (jii > 2.65 Mp. 

Now we consider this situation for the a-model. From 
Fig. 2 we see the fine tuning can be ameliorated con¬ 
siderably. For n = 0 (n = 3) the minimum field value 
to achieve = 0.7 and wq < —0.9 is (j)i > 1.95 Mp 
(0.9 Mp) for a = 1. The conditions to reach wq closer to 
— 1 are even more softened, since once the a-model field 
is on the potential plateau then w can easily be extremely 
close to —1. Again, the plateau effectively stretches out 
hilltop potentials, such as from pseudo-Nambu goldstone 
boson or axion fields, or Higgs symmetry breaking, vastly 
increasing the initial conditions delivering w ~ — 1. 

We can now further generalize the a-model by taking 

(where before we restricted to p = 2). This now acts 
at small (j) like any chosen monomial potential (p^, while 
retaining the plateau at larger field values. True attrac¬ 
tors to ic = — 1 can now be achieved when the field is 
at values x > p/{2n — p), requiring n > p. For exam¬ 
ple, rather than matching the quadratic potential near 
the minimum we can use the axion monodromy poten¬ 
tial i.e. p = 2/3. With n = 1, say, this has a late 
time de Sitter attractor for x > 0.5 or pi > 1.35 i/a. 

Let us explore the attractor behavior in the inflation¬ 
ary sense of an a-attractor, where there is a common be¬ 
havior for observables for some limit of a. When a ^ 1, 
all models regardless of n will look like the V oi pP case, 
e.g. a quadratic potential for our baseline of p = 2. Con¬ 
versely, when a —^ 0, all models will see a plateau poten¬ 
tial. In particular, when n > p and the held is beyond the 
potential maximum at a; = p/(2n — p), then all models 
are attracted to the de Sitter state. 

Figure 4 shows the dynamics in the w — w' phase space, 
where w' = dw/dlna, for three values of p, hxing n = 3 
and taking pi = 1.5-^a. For a = 1, the dynamics are 
well separated, though with common general character¬ 
istics. Note that they all start off as thawing dark en¬ 
ergy, evolving along the canonical thawing behavior of 
w' = 3(1-I-w), and even at present (denoted by x’s) they 
lie in the thawing region, at roughly w' = l-|-w. However 
in the future they turn around and convert to freezing be¬ 
havior, heading along an attractor toward the de Sitter 
point w = — 1, u>' = 0. For small a, these characteris¬ 
tics persist but curves for all p stay much closer to the 
ultimate cosmological constant behavior. 

Returning to the baseline (p = 2) case, if we expand 
the a-model potential for small /)/i/a, we see 

V{p^y/a)Ki ^p^ (^-^=p + 0{p^/a:^ . (13) 

So there is an attractor behavior for large a where for 
any n the potential acts like a simple quadratic potential. 
Figure 5 illustrates this behavior. (An analogous result 
will hold for any p.) 



FIG. 4. The dynamics in the w — w' phase space transitions 
between thawing and freezing dark energy, with an attrac¬ 
tor behavior both toward a cosmological constant and toward 
common dynamics as a —> 0. Models with p = 2/3, 1, and 
2 (red, blue, magenta, or from outer to inner) are shown, for 
0 = 1 (solid curves) and a = 0.1 (dotted curves). The present 
is denoted by x’s and the curves are stopped at a = 220 for 
clarity. Note the p = 2, a = 0.1 case has not quite turned 
around since it is very close to the potential maximum. 


As an interesting additional point, Eq. (13) shows that 
away from the minimum the potential becomes anhar- 
monic. That is, the potential slope becomes shallower, 
just as for hilltop potentials. We know these have neg¬ 
ative mass squared terms near the maximum (also see 
[9-11] for other anharmonic instabilities), and we find 
this as well for the a-model. Specifically, 

= 1/" = ^ (TT^ ^ 

—x^{Z — n — 2n^) — x^{Z — 5n + 2'n?)\ . 

Note that the mass vanishes for 1, as do all the 

derivatives of the potential; the plateau becomes flat. 
Figure 6 exhibits the mass squared as a function of p 
for various n. 

For hilltop potentials the spinodal, or tachyonic, insta¬ 
bility when the mass squared is negative is not a major 
concern [12-14], but for a-models it does exist for a wider 
range of field values. This raises the possibility of inter¬ 
esting phenomenology, such as clustering dark energy. 
Here we consider only the background, zero momentum 
modes, but will investigate this further in future work. 








5 



W 


FIG. 5. The dynamics in the w — w' phase space increasingly 
resembles a quadratic potential (dashed, long magenta curve) 
as a —>■ 00 . Models with n — 0, 1, and 2 (black, blue, red, or 
from longest to shortest curves for each line type) are shown, 
for a = 1 (solid curves) and a = 10 (dotted curves). For a = 
10, the curves, and the endpoints at a = 1, move noticeably 
closer to the quadratic result. 
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FIG. 6. The second derivative of the potential, giving the 
effective mass squared, is plotted as a function of the field 
(j) in units of ^/a, for various values n (increasing from top 
to bottom curves at small cj)). Note that all derivatives ap¬ 
proach zero at large field values (since the potential plateau 
is asymptotically flat). 


V. CONCLUSIONS 


We investigated the properties of a-attractors, as re¬ 
cently highlighted for inflation, as dark energy for late 
time acceleration instead. Such models have several in¬ 
teresting properties, combining aspects of the thawing 
and freezing classes of dark energy together. By general¬ 
izing the usual a-attractors we derive models that have 
a true de Sitter attractor, as well as ones that have a 
metastable acceleration with equation of state ratio lin¬ 
gering near w « — 1. 

The dynamics near the potential minimum, where the 
Qf-model looks like a quadratic potential, is found to be 
less fine tuned that a standard quadratic potential, espe¬ 
cially for values of w consistent with observations. Fur¬ 
ther from the minimum, the model resembles a stretched 
hilltop model, with elements of an exponential potential 
as well. 

A family of theories can be defined as one varies the 
value of a. As a —>■ 0, the field sees predominantly a 
plateau, with a slight slope up or down, depending on the 
value of the generalized parameter n. This determines 
whether there is a true de Sitter attractor or not, but 
either way w k, for many e-folds of expansion. For 
a ^ 1, the potential increasingly resembles a quadratic 
potential, or any monomial (jf in a further generalization. 

We note the a-model has anharmonic properties, and 
will have negative effective mass squared in some regions. 
Either of these can induce clustering in the scalar field, 
possibly leading to interesting effects. While here we con¬ 
centrate on the background field, future work will explore 
the possibility of observable signatures of these fluctua¬ 
tion properties. 
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